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SECTION  I 


INTRODUCTION 


Maneuverable  supersonic  and  hypersonic  flight  vehicles 
require  a  variety  of  control  surfaces  that  generate  complex 

flow  patterns.'  One  common  example  of  such  a  control  surface 
is  a  blunt  fin,  which  at  locally  supersonic  flow  conditions  can 
induce  a  shock/boundary  layer  interaction  associated  with  flow 

separation.  The  baseline  case  for  the  blunt  fin  control 

surface  is  the  blunt  fin  on  a  flat  plate  at  zero  angle  of 
attack  to  the  freestream,  as  shown  in  Figure  1.  Although  this 
configuration  has  been  previously  investigated  both 
experimentally  and  numerically,  much  remains  unknown  about  the 
flow  field,  especially  at  non-zero  sweepback  angles.  The 
objective  of  the  present  research  is  to  establish  computational 
evidence  describing  the  effect  of  the  sweepback  angle.  A,  on 
the  number  and  strength  of  the  upstream  vortical  structures 
associated  with  diverting  flow  past  the  fin.  For  the  case  of 

the  vertical  fin  (A=0) ,  Hung  et  al  '  predicted  the  nature  of 
the  flow  field  structure  using  numerical  studies  based  on 

previous  experiments.^  ' Their  results  showed  the  formation  of 
a  pair  of  horseshoe  vortices  upstream  of  the  fin,  with  the 
primary  vortex  diverting  entrained  freestream  fluid  past  the 
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fin.  Other  numerical  work  by  Shang  and  Scherr,  however, 
indicated  that  a  highly  swept  fin  on  a  complex  configuration 
had  no  associated  upstream  horseshoe  vortex  to  divert  the 

fluid.  Earlier  experimental  evidence  ’  also  indicated  that 
sweepback  had  an  effect  in  reducing  the  upstream  influence  of  a 
fin,  although  no  firm  conclusions  were  made  about  the  existence 
or  strength  of  vortical  structures. 

The  above  results  indicate  that  the  flow  structure  should 
undergo  a  bifurcation  at  a  sweepback  angle  sufficiently  high 
that  the  horseshoe  vortex  upstream  of  the  fin  leading  edge 
ceases  to  exist.  At  this  juncture,  an  associated  change  must 
take  place  in  the  topological  phase  portrait  on  the  surface  of 

the  flat  plate. ^  In  addition,  the  decrease  in  upstream 
influence  associated  with  increased  fin  sweepback  should  also 
coincide  with  a  less  pronounced  spanwise  propagation  of  flow 
disturbances.  This  information  is  important,  since  vortical 
structures  can  cause  significant  variations  in  pressure  and 
rate  of  heat  transfer  on  neighboring  surfaces  as  well  as  on  the 
blunt  fin  itself,  thus  altering  control  effectiveness  in  the 
entire  region  near  the  fin. 

A  two  step  procedure  was  followed  to  determine  the  effect 
of  fin  sweepback  on  the  upstream  topological  structure  and 
spanwise  propagation  of  flow  disturbances.  Results  were  first 
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calculat^ed  echoing  the  vertical  blunt  fin  flow  conditions  used 

by  previous  numerical  and  experimental  investigators.^  ^  The 

stagnation  temperature  and  pressure  were  260  K  and  6.8  x  10^ 

Nm  ,  respectively.  The  freestreaun  Mach  number  was  2.95, 
yielding  a  unit  Reynolds  number  of  64  million  per  meter.  The 
boundary  layer  was  fully  turbulent  with  a  thickness  of  1.27  cm 
at  the  fin  leading  edge.  Th t  fin  loading  edge  diameter  was 
also  1-27  cm.  The  compu tat. i ons  were  accomplished  using  a 
tlirec  dimensional  finite  difference  explicit  code  based  on 

MacCormac.k  ’  s  1969  predictor-corrector  algorithm.'^  Once  the 
code  was  verified  for  the  vertical  fin,  it  was  utilized  under 
the  same  upstream  flow  conditions  for  increasing  fin  sweep  to 
assess  the  formation  and  strength  of  the  vortical  structures  as 
a  function  of  the  sweep  angle.  This  procedure  was  continued 
until  the  bifurcation  point,  reflected  by  a  change  in  the  phase 
port-rait  of  the  surface  shear  stress  vector,  was  reached. 


SECTION  II 


ANAI.YSIS 
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1 .  TOPOLOGICAL  BIFURCATION 


Any  'three-dimensional  flow  field  must  satisfy  topological 
rules  which  arise  from  the  requirement  of  flow  continuity.  A 
separated  flow  can  be  characterized  by  the  number  of  node  and 
saddle  points  on  a  given  plane  cutting  tlie  tliree  dimensional 
flow  field;  traditionally,  the  flow  adjacent  to  a  surface  is 
represented  by  the  skin- f ri ction  line  pattern  formed  by 
experimental  oil  streaks.  The  status  of  the  surface  shear  flow 
pattern  is  termed  the  phase  portrait  of  the  surface  shear 
stress  vector  (path -preserv i ng  mapping) .  If  the  character  of 
the  phase  portrait  changes  due  to  a  charige  in  some  external 
flow  parameter,  the  phase  j)cji-trait  is  said  to  undergo 

bifurcation.^  The  flow  past  a  vertical  blunt  fin  is 
characterized  by  a  surface  oil  pattern  on  the  flat  plate  that 
includes  a  convergence  line  (indicating  that  flow  leaves  the 
surface  to  roll  up  into  the  horseshoe  vortex)  and  a  divergence 
line  (indicating  reattachment)  upstream  of  the  fin  leading 
edge.  This  pattern  is  seen  in  the  numerically  simulated  oil 
patterns  of  Figure  2,  which  was  taken  from  Reference  7. 


In  contrast  to  the  pattern  seen  for  the  vertical  fin. 
Figure  3  shows  the  surface  shear  pattern  near  a  highly  swept 
fin  on  a  complex  body,  as  computed  in  Reference  8.  Here,  the 
absence  of  a  set  of  convergence  and  divergence  lines  upstream 


of  the  fin  indicates  fully  attached,  vortex-free  flow.  Thus 
given  a  flow  represented  by  the  surface  shear  portrait  shown  in 
Figure  2,  if  one  were  to  increase  the  fin  sweepback  angle,  at 
some  sufficiently  large  sweepback  angle  the  surface  shear 
pattern  should  bifurcate.  The  surface  shear  portrait  would 
then  be  similar  to  that  in  Figure  3,  characterized  by  the 
absence  of  distinct  convergence  and  divergence  lines  upstream 
of  the  fin. 


2.  GOVERNING  EQUATIONS 


The  flow  under  consideration  can  be  described  by  the  fully 
three  dimensional,  mass  averaged  Navier-Stokes  equations  with 
an  algebraic  turbulence  model: 
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7he  assumption  is  made  that  the  mass  averaged  equations 
represent  the  ensemble  average  of  the  flow,  even  though  the 

flow  is  unsteady  due  to  shock  oscillation.^  7he  frequency  of 
the  shock  oscillation  is  in  the  kHz  range,  and  the  assumption 
is  supported  by  the  referenced  experimental  findings  which 
exhibit  no  unsteadiness  in  the  measurements  of  the  mean  flow 
properties . 
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3 .  NUMERICAL  PROCEDURES 


The  numerical  scheme  for  the  present  study  is  MacCormack’s 
1969  predictor-corrector  method,  which  has  been  shown  to  yield 
accurate  results  for  high-Reynolds  number  supersonic  flows. 

The  computational  grids  were  algebraically  generated,  and 
were  dimensioned  40  x  32  x  32  in  the  streamwise,  spanwise,  and 
vertical  directions,  respectively.  The  grid  structure  for  the 
vertical  fin  (A=0*)  is  very  similar  to  that  used  in  References 
6  and  7.  However,  the  minimum  grid  spacing  of  .00433D  used  in 
the  referenced  numerical  work  was  found  to  be  too  coarse  for 
the  present  purpose,  because  it  yields  a  value  of  y+  at  the 
wall  of  ~30.  While  this  resolution  may  be  sufficient  to 
capture  the  major  structures  of  the  flow  past  the  vertical  fin, 
since  those  structures  are  predominately  inviscid,  the  ability 
to  predict  the  upstream  flow  field  near  the  fin-flat  plate 
juncture  for  configurations  of  high  sweepback  was  seriously 
degraded  by  inadequate  resolution  of  the  viscous  region. 
Therefore  a  minimum  grid  spacing  of  .OOlD  was  used  at  the  flat 
plate  surface  for  all  the  grids  that  were  dimensioned  40  x  32  x 
32.  This  resulted  in  a  y+  at  the  wall  of  ~10.  In  the 
direction  i.ormal  to  the  fin,  however,  the  minimum  spacing  at 
the  fin  surface  was  kept  at  . 004D .  This  created  a  less  severe 
stretching  in  the  normal  direction,  which  maintained  sufficient 
resolution  upstream  of  the  fin  leading  edge  near  the 
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experimentally  determined  location  of  the  convergence  of  the 
surface  shear-stress  lines.  Figure  4a  is  an  elevated  view  of 
the  grid  used  to  compute  the  flow  past  the  vertical  fin.  The 
sweptback  grids  were  generated  using  the  grid  for  the  vertical 
fin  by  moving  each  horizontal  plane  rearward  an  appropriate 
amount  according  to  the  sweepback  angle.  Figure  4b  shows  the 
elevated  view  of  the  grid  for  the  60*  swept  fin.  The  severity 
of  the  stretching  of  individual  cells  for  this  grid  can  be  seen 
in  Figure  4c,  which  is  a  closeup  of  the  grid  cells  near  the 
fin-flat  plate  junction  on  the  plane  of  symmetry. 

The  turbulence  model  was  the  Baldwin-Lomax  eddy  viscosity 
model,  with  a  modified  scaling  length  designed  to  take  into 
account  the  presence  of  two  walls.  In  this  formulation,  the 
scaling  length  is 
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where  z  and  y  are  the  normal  distances  from  the  two  walls.  The 
general  formulation  of  the  eddy  viscosity  near  two  walls  was 

considered  in  detail  in  previous  references^*'' and  is  not 
explained  here. 
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The  boundary  conditions  were  essentially  the  sajne  as  those 
used  in  Reference  7.  The  upstream  boundary  condition  was 
prescribed  as  the  fully  turbulent  flat  plate  boundary  layer  of 
thickness  1.27  cm,  matching  the  fin  leading  edge  diameter. 

This  was  obtained  by  computing  the  flow  on  a  flat  plate  until 
convergence,  and  using  the  station  with  the  appropriate 
boundary  layer  thickness  as  the  upstream  plane  for  calculation 
of  the  blunt  fin  flow.  The  surface  boundary  conditions  were 
the  standard  no-slip,  adiabatic  wall  conditions  with  zero 
normal  pressure  gradient.  Zero  gradient  conditions  were 
applied  on  the  downstream  and  normal  boundaries,  to  ensure  the 
flow  remained  unperturbed  along  the  fin  and  at  large  distances 
from  the  intersection  of  the  surfaces. 


SECTION  III 

RESULTS 

1.  FLAT  PLATE  OIL  STREAKS 

Simulation  of  oil  streaks  on  the  flat  plate  surface  is 
one  method  of  evaluating  the  upstream  and  spanwise  effects  of 
the  presence  of  the  blunt  fin  in  the  flow  field.  Figure  5 
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presents  the  computed  streak  patterns  from  the  present  study 
for  the  case  of  the  vertical  fin.  The  reference  lines  on  the 
flat  plate  are  1/2  diameter  apart.  The  similarity  between 
these  results  and  those  shown  in  Figure  2  is  clear.  The 
primary  convergence  line  for  the  present  calculation  is 
slightly  more  than  two  diameters  upstream  of  the  fin  leading 
edge,  and  does  not  recede  downstream  as  quickly  as  the  computed 
convergence  line  shown  in  Figure  2.  In  this  respect  the 
current  calculation  more  closely  resembles  the  experimental 
results,  shown  as  the  dashed  line  in  Figure  2.  The  extent  of 
spanwise  disturbance  can  be  judged  from  the  location  of  the 
convergence  line  as  it  recedes  downstream  and  outward  past  the 
fin.  At  the  fin  leading  edge,  the  spanwise  location  of  the 
limiting  streamline  is  3.8  diameters,  increasing  to  5.2  and 
then  6.3  diameters  at  stations  two  and  four  diameters 
downstream  of  the  leading  edge,  respectively.  Figure  6  shows 
similar  views  of  the  oil  patterns  computed  for  cases  of  fin 
sweepback  of  30* ,  45* ,  60* ,  and  68* .  The  upstream  location  of 
the  convergence  line  approaches  the  fin  steadily  as  the 
sweepback  increases,  although  there  is  a  distinct  accumulation 
line  upstream  of  the  fin  even  at  60*. 

At  the  sweepback  angle  of  68*,  however,  no  convergence  line 
can  be  found  upstream  of  the  leading  edge.  This  corresponds  to 
the  bifurcation  point:  a  new  phase  portrait  of  the  surface 
shear-stress  vector  has  emerged.  There  is  a  simultaneous 
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change  in  the  external  flow  pattern  associated  with  the  phase 
portrait.  The  horseshoe  vortex,  which  previously  was  the 
mechanism  for  diverting  oncoming  fluid  around  the  fin,  no 
longer  exists.  What  may  appear  to  be  a  convergence  line 
spreading  spanwise  from  the  leading  edge  of  the  fin  is  actually 
a  cluster  of  parallel  streamlines  which  are  forced  into  close 
proximity  due  to  the  deflection  of  the  flow  originating 
upstream  of  the  fin.  The  region  between  these  compressed 
streamlines  and  the  blunt  fin  is  occupied  by  fluid  which  passes 
over  the  top  of  the  fin  before  sweeping  down.  By  comparing  the 
location  of  the  computed  oil  patterns  with  those  for  the 
previous  sweepback  angles,  one  can  see  that  the  extent  of  the 
disturbed  flow  is  very  small  for  this  case.  This  decreasing 
spanwise  propagation  of  disturbance  is  in  a  consistent  trend 
with  the  assertion  that  increasing  fin  sweepback  indicates 
weakening  of  the  horseshoe  vortex.  The  effect  of  the  sweepback 
on  the  spanwise  location  of  the  oil  accumulation  line  is 
summarized  in  Figure  7,  which  plots  the  convergence  line 
location  versus  sweep  angle  for  the  range  of  zero  to  sixty 
degree  sweep. 

2.  FLAT  PLATE  PRESSURE  DISTRIBUTION 

Although  the  location  of  the  oil  convergence  line  gives 
an  indication  of  where  the  flow  leaves  the  surface  to  roll  up 
into  vortical  structures,  it  does  not  specifically  define  the 
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extent  of  fin  influence  on  the  flow  field.  One  other  parameter 
which  can  help  define  this  influence  is  the  distribution  of 
pressure  on  the  flat  plate.  The  surface  pressure  distribution 
in  the  plane  of  symmetry  for  the  vertical  fin  is  presented  in 
Figure  8.  This  pressure  distribution  reveals  the  strong 
influence  exerted  by  the  embedded  horseshoe  vortical  structure. 
Twin  pressure  peaks  are  partitioned  by  a  trough  which  is 
induced  by  the  presence  of  a  strong  primary  vortex  system 
attached  to  the  surface.  The  present  calculation  produces  good 
agreement  with  previous  experiments  and  calculations.  The  dip 
in  pressure  near  the  leading  edge  of  the  blunt  fin  is  similar 
to  the  dip  observed  in  the  calculations  of  Reference  7  when  a 
refined  grid  was  used. 

The  particular  behavior  of  the  surface  pressure  at  the 
plane  of  symmetry  persists  beyond  the  blunt  leading  edge  of  the 
fin.  The  twin  peaks  characterizing  the  distribution  are 
diminished  by  the  rapid  expansion  around  the  leading  edge  and 
by  the  dissipation  of  the  horseshoe  vortex  as  the  flow  moves 
downstream.  This  process  can  be  observed  in  the  contours  of 
pressure  on  the  flat  plate,  shown  in  Figure  9  for  both 
experimental  and  computational  results.  Examination  of  Figures 
9a  and  b  reveal  an  excellent  agreement  between  experiment  and 
computation  for  the  vertical  fin.  The  initial  limit  of 
disturbance  is  approximately  2.7  diameters  upstream  of  the  fin, 
receding  to  less  than  one  diameter  upstream  of  the  fin  at  a 
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spanwise  station  of  four  diameters.  The  first  peak  pressure 
upstream  of  the  fin  leading  edge  is  in  the  vicinity  of  1.3 
diameters  upstream,  although  the  peak  is  somewhat  smeared  in 
the  calculation.  The  location  of  the  pressure  trough 
downstream  of  the  first  peak  corresponds  to  the  experimentally 
determined  location.  The  propagation  of  the  trough  and  the 
expansion  around  the  fin  also  compare  very  well  with 
experimental  measurements.  In  Figure  9c,  A=30* ,  one  can 
observe  that  the  initial  pressure  disturbance  has  receded  to 
one  diameter  upstream  of  the  fin  leading  edge.  The  first  peak 
pressure  retains  the  same  magnitude  as  that  for  the  vertical 
fin,  but  is  located  much  closer  to  the  fin.  The  pressure 
trough  also  doesn’t  change  significantly  in  magnitude,  although 
it  does  change  in  location.  The  final  peak  pressure  on  the  fin 
is  lower  than  the  counterpart  on  the  vertical  fin.  The 
spanwise  propagation  of  flow  disturbance  is  also  lessened 
significantly,  as  the  pressure  disturbance  does  not  reach  a 
spanwise  station  of  four  diameters  until  more  than  two 
diameters  downstream  of  the  fin  leading  edge.  For  cases  of 
sweepback  of  45,  60,  and  68  degrees,  the  first  peak  pressure 
upstreaiD  of  the  fin  leading  edge  remains  remarkably  constant. 
The  pressure  trough  between  the  first  and  final  pressure  peaks 
becomes  less  pronounced  with  greater  sweepback,  however,  while 
the  final  peak  pressure  decreases  in  magnitude.  The  result  is 
that  at  a  sweepback  between  60  and  68  degrees,  the  pressure 
trough  disappears  entirely  and  there  is  only  one  pressure  peak. 
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This  coincides  with  the  bifurcation  of  the  surface  shear 
portrait  and  the  disappearance  of  the  horseshoe  vortex  upstream 
of  the  fin.  The  spanwise  disturbance  also  continuously 
decreases,  until  at  60  degrees,  the  disturbance  reaches  four 
spanwise  diameters  at  a  downstream  station  of  five  diameters. 

3 .  VORTICAL  STRUCTURES 

The  primary  mechanism  associated  with  diverting  fluid 
past  the  blunt  fin  is  a  horseshoe  vortex  which  forms  upstream 
of  the  f  in.  The  structure  of  this  vortex  for  the  vertical  fin 
in  the  plane  of  symmetry  can  be  seen  by  the  particle  traces  of 
Figure  10.  In  this  and  following  figures,  the  reference 
coordinate  lines  on  the  flat  plate  are  1/2  diameter  apart.  The 
center  of  the  vortex  is  located  slightly  more  than  1/2  diameter 
upstream.  This  corresponds  to  the  location  of  the  pressure 
trough  upstream  of  the  blunt  fin  in  Figure  9b.  Figure  11  shows 
the  overhead  view  of  the  trajectory  of  particles  swept  into 
this  vortex.  Again  comparison  with  Figure  9b  shows  a  strong 
correlation  between  the  vortex  trajectory  and  the  pressure 
trough  present  in  the  flow  field.  For  example,  at  a  distance 
of  two  diameters  downstream  of  the  fin  leading  edge,  the 
computed  vortex  particle  traces  pass  a  location  at  three 
diameters  outboard  of  the  center  line.  This  coincides  with  the 
experimental  and  computed  location  of  the  pressure  trough . 
Figures  12-14  show  further  evidence  to  substantiate  this 
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correspondence.  Note  that  in  these  figures,  the  blunt  fin  is 
not  shown  for  clarity.  Figure  12  shows  the  vortex  trajectory 
for  A=30*  .  At  a  station  1/2  diaxoeter  downstresim  of  the  fin 
leading  edge,  the  center  of  the  vortex  is  located  at  a  distance 
of  one  diameter  from  the  center  line.  This  coincides  closely 
with  the  location  of  the  pressure  defect  as  plotted  in  Figure 
9c.  The  center  of  the  vortex  in  the  symmetry  plane  is  located 
at  roughly  .2  diameters  of  the  leading  edge  of  the  fin,  which 
again  closely  correlates  to  the  location  of  the  pressure  trough 
in  Figure  9c.  From  Figure  13,  A=45* ,  one  can  see  the  vortex 
situated  a  very  small  distance  upstream  of  the  fin,  and  at  one 
diameter  downstream,  the  vortex  is  located  slightly  more  than 
one  diameter  spanwise  from  the  center  line.  These  results 
again  agree  with  the  computed  location  of  the  pressure  trough 
(Fig.  9d) .  As  the  sweepback  angle  becomes  greater,  the  vortex 
gradually  decreases  in  size,  which  also  corresponds  to  the 
weakening  of  the  pressure  trough  shown  with  increasing 
sweepback  angle.  Figure  14  shows  the  continuation  of  this 
trend.  Here,  the  vortex  forms  directly  above  the  fin  leading 
edge  and  is  still  located  within  one  spanwise  diameter  of  the 
blunt  fin  when  it  is  one  diameter  downstream  of  the  leading 
edge.  By  this  time  the  dimension  of  the  vortex  is  smaller  than 
l/20th  the  boundary  layer  thickness  in  diameter  and  shows 
little  effect  on  the  mean  flow.  This  observation  can  be 
confirmed  by  the  pressure  contour  plot  Fig.  9d ,  which  shows 
only  a  slight  wiggle  in  the  contours  near  the  intersection  of 
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the  fin  and  the  plate  to  indicate  the  existence  of  the  pressure 
drop.  Figure  15  is  an  elevated  view  of  the  flow  field  very 
near  the  fin-flat  plate  juncture  for  A=68* .  It  is  apparent,  as 
stated  before,  that  no  vortex  formation  is  present.  Upstream 
fluid  particles  approaching  the  fin  simply  divert  to  the  side 
of  the  fin  without  being  entrained  by  a  vortical  structure,  or 
continue  up  the  fin  before  passing  down  its  side. 


SECTION  IV 

CONCLUSION 


The  flow  past  a  blunt  fin  at  Mach  2.95  was  successfully 
calculated  for  several  fin  sweepback  angles.  Numerical  results 
show  that  the  flow  field  topological  structure  undergoes  a 
bifurcation  at  sufficiently  high  angles  of  sweepback. 

Sweepback  of  the  fin  reduces  the  upstream  influence  of  the  flow 
field  and  simultaneously  weakens  the  horseshoe  vortex  upstream 
of  the  fin.  Spanwise  influence  is  thereby  decreased.  For  the 
simulated  conditions,  the  horseshoe  vortex  ceased  to  exist 
between  60  and  68*  sweep . 
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Figure  I 


The  Baseline  case  for  the  Blunt  Fin  Control  Surface  is  the  Blunt 
Fin  on  a  Flat  Plate  at  Zero  Angle  of  Attack  to  the  Freestream 
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Figure  5:  The  Computed  Streak  Patterns  From  the  Present  Study  for  the  Case  of  The  Vertical  Fin 


Figure  7:  Plots  the  Convergence  Line  Location  Versus  Sweep  Angle  for  the 
Range  of  Zero  to  Sixty  Degree  Sweep 
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Preseat  Compatatleas 
Compatstioas  (Ref  7) 
Experimeat  (Ref  4) 


The  Surface  Pressure  Distribution  in  the  Plane  of  Symmetry  for  the 
Vertical  Fin 
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erting  Fluid  Past  the  Blunt  Fin  is  a  Horseshoe  Vortex  Which  Forms  Upstream 
the  Fin 
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Figu 
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Figure  14:  Shows  the  Continuation  of  the  Trend 


34 


Figure  15:  Elevated  View  of  the  Flowfield  Very  Near  the  Fin-Flat  Plate  June 


